Substrate-induced antiferromagnetism of an Fe monolayer on the Ir(OOl) surface 
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We present detailed ab initio study of structural and magnetic stability of a Fe-monolayer on the 
fcc(OOl) surface of iridium. The Fe-monolayer has a strong tendency to order antiferromagnetically 
for the true relaxed geometry. On the contrary an unrelaxed Fe/Ir(001) sample has a ferromagnetic 
ground state. The antiferromagnetism is thus stabilized by the decreased Fe-Ir layer spacing in 
striking contrast to the recently experimentally observed antiferromagnetism of the Fe/W(001) 
system which exists also for an ideal bulk-truncated, unrelaxed geometry. The calculated layer 
relaxations for Fe/Ir(001) agree reasonably well with recent experimental LEED data. The present 
study centers around the evaluation of pair exchange interactions between Fe-atoms in the Fe- 
overlayer as a function of the Fe/Ir interlayer distance which allows for a detailed understanding of 
the antiferromagnetism of a Fe/Ir(001) overlayer. Furthermore, our calculations indicate that the 
nature of the true ground state could be more complex and display a spin spiral-like rather than a 
c(2x2)-antiferromagnetic order. Finally, the magnetic stability of the Fe monolayer on the Ir(OOl) 
surface is compared to the closely related Fe/Rh(001) system. 

PACS numbers: 75.30.-m,81.20.-n 

Keywords: surface magnetism, iron, iridium, density functional calculation 



I. INTRODUCTION 



Magnetic overlayers, i.e., the thin films of magnetic 
materials on the nonmagnetic substrate are systems 
with great technological potential. Magnetic overlay- 
ers are also a convenient system for a deeper under- 
standing of the origin of magnetism in the solid state. 
The prototypical system is a magnetic monolayer (ML) 
on a non-magnetic substrate which was the subject 
of many theoretical and experimental studies in the 



past 



1,2,3,4,5,6,7,8,9,10,11.12,13 



Yet, a technological prepara- 



tion, experimental study and detailed theoretical under- 
standing of the overlayer magnetism is still a challenge 
to solid state physics. 

There are a few important features which distinguish 
magnetic overlayers from conventional magnetic materi- 
als. This is, first of all, the presence of an external agent, 
namely the presence of a non-magnetic substrate, which 
can strongly modify the magnetic state and properties 
of the overlayer system: (i) The reduced coordination on 
the surface induces geometry as well as chemical binding 
changes, (ii) Magnetic atoms adapt to the underlying 
lateral substrate lattice spacings, particularly true for a 
monolayer case, (iii) The substrate electronic structure, 
in particular the position of the substrate Fermi level is 
relevant for overlayer magnetism, (iv) The position of 
Fermi level can be tuned by the substrate electron con- 
centration, e.g. by growing the overlayer on an alloy 
substrate with varying concentration of the constituents. 
Alternatively, one could consider a substrate alloy which 
contains magnetic and non-magnetic species or only mag- 
netic species, and/or a partial coverage of the substrate 
to vary the properties of the magnetic overlayer. Theo- 
retical and experimental studies on the above trends help 



to establish a deeper understanding of the origin of sur- 
face magnetism. 

First-principles calculations represent a powerful tool 
for such studies, as they allow to determine reliably the 
underlying lattice structure (possible layer relaxations, 
surface reconstructions, etc.) which can be directly com- 
pared with experiment (LEED). They also allow to find 
out the underlying magnetic structure although in this 
respect the situation is much more complex. In addi- 
tion to collinear magnetic configurations (ferromagnetic 
(FM) and antiferromagnetic (AFM) ones) more complex 
configurations may exist (e.g. recently observed chi- 
ral structures in bcc-Fe/W(001) 1 , magnetism of random 
overlayers 2 -, etc.). Consequently, an estimation and dis- 
cussion of exchange interactions is very useful to gain a 
deeper understanding of properties of both overlayer and 
bulk magnets and magnetic alloys, including the diluted 
magnetic semiconductor s 14 ' 15 . Surprisingly and in con- 
trast to bulk systems, studies of exchange interactions for 
magnetic overlayers are still very rare 3 ^ despite of their 
obvious importance. In particular, the distance depen- 
dence of exchange integrals can be very different from 
that in the bulk if, for example, adatoms can interact via 
the host surface state as e.g. Co-adatoms on the fcc(lll) 
faces on noble metals 4 ^. 

Great emphasis, both theoretical and experimental, 
has been put recently on the study of the magnetic prop- 
erties of Fe-overlayers on the (001) and (110) faces of 
bec tungsten. An unusual AFM state was predicted 
theoretically^ and confirmed experimentally for the bec- 
Fe/W(001) system^ (in fact, the true ground state seems 
to be a more complex one exhibiting a chiral stated). 
There is also evidence from theory that the ground state 
of the Co/W(001) overlayer is AFM, whereas Mn and Cr 
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overlayers have a FM ground stated. This means that a 
(Fe,Mn)-alloy overlayer on the bcc-W(OOl) should exhibit 
a crossover from an AFM to a FM ground stated. Fi- 
nally, a similar AFM- to FM-crossover was predicted for 
the bcc-Fe/(Ta,W)(001) alloy substrate system^ ., where 
the ground state of the bcc-Fe/Ta(001) is FM. However, 
only the leading exchange integrals were estimated in the 
latter case^. 

The (OOl)-faces of the fee- and bec-substrates exhibit a 
simple rectangular array of lattice sites and differ essen- 
tially by the d : a^-ratio, where is the lateral lattice 
constant and d is the interlayer distance. Therefore an 
investigation of the magnetic properties of Fe-overlayers 
on the (001) faces of fee transition metals poses an in- 
teresting problem. A particularly suitable system is the 
fcc-Fe/Ir(001) where very thin overlayers were grown suc- 
cessfully with negligible Fe-Ir intermixing 11 . Further- 
more, there exist reliable LEED measurements^ eluci- 
dating the detailed geometry, while preliminary MOKE- 
studies 11 indicate no magnetic signal in the limit of 
a monolayer coverage; a situation similar to the bec- 
Fe/W(001). 

The clean Ir(001) surface undergoes the (5 x 1) quasi- 
hexagonal reconstruction. The finite Fe-coverage (larger 
than 0.25 monolayer) lifts the reconstruction 1 -. In order 
to avoid possible corrugation and Fe-Ir intermixing the 
metastable unreconstructed (1 x l)-Ir(OOl) surface was 
prepared which is, however, stable for room temperatures 
and below (for details see Ref . fill) . 

The aim of the present study is a comprehensive 
first-principles investigation of the properties of an fec- 
Fe/Ir(001) monolayer system. In a first step we per- 
form a structural study of Fe/Ir(001) by using two 
highly accurate DFT methods, namely the WIEN2kA& 
and VAS P 17 i 18 codes. In the next step we investigate the 
magnetic structure of the geometrically relaxed system 
by comparing total energies of the non-magnetic (NM) , 
FM-, and AFM-configurations (c(2x2)-AFM). Here, the 
main problem is the large number of possible magnetic 
configurations (including the non-collinear ones) to be 
considered. Instead of a brute-force search, we evaluate 
the total energy of the disordered local moment (DLM) 
state similarly to Ref.[T3- The DLM state is a model state 
with zero total magnetic moment which results from the 
disorder of spin orientations of otherwise non-zero local 
magnetic moments^ A strong indicator for a more com- 
plex magnetic state of the system under consideration is 
a lower total energy of the DLM state as compared to the 
NM and FM states. One should note, however, that if 
some specific magnetic state is the ground state, e.g., the 
c(2x 2)-AFM, its total energy is usually lower than that of 
the DLM state. On the other hand, in random magnetic 
systems, non-collinear states can be the ground state 
as found, e.g. in fcc-NiMn alloys22 or in (Cu,Ni)MnSb 
Heusler alloys—. The usefulness of the DLM concept was 
demonstrated recently for both overlayer studies^ and 
disordered magnetic semiconductors^. The DLM pic- 
ture may be straightforwardly implemented in the frame- 



work of the coherent-potential approximation (CPA) 22 . 
Therefore in a third step, we perform studies based on 
the Green function implementation of the TB-LMTO- 
CPA method 23 in the framework of the surface Green 
function (SGF) approach, in addition to the above men- 
tioned collinear WIEN2k and VASP calculations. The 
TB-LMTO-SGF approach employs a realistic semiinfi- 
nite sample geometry (no slabs or periodic supercells) 
and allows to implement the DLM model. The one- 
electron potentials are treated within the atomic sphere 
approximation (ASA); the dipole barrier due to the sam- 
ple electrons in the vacuum is included in the formal- 
ism. TB-LMTO-SGF even allows to include the effect 
of layer relaxations 10 , provided, they are known either 
from full-potential calculations or from experiment. An 
important advantage of the TB-LMTO-SGF approach is 
the possibility to estimate exchange interactions between 
magnetic atoms in the overlayer by a straightforward gen- 
eralization of the well-approved bulk concep t 15 ' . Sum- 
marizing, the TB-LMTO-SGF approach is a very useful 
tool for a qualitative understanding of the results while 
the full-potential approaches are superior concerning the 
quantitative values and thus provide so to say the corner- 
stones. We will demonstrate that at least in the present 
case the results of both types of calculations are in a good 
quantitative agreement which justifies our assumptions 
and reinforces our conclusions. 



II. COMPUTATIONAL DETAILS 

First principle density functional theory calculations 
were performed using both the all-electron full poten- 
tial linearized augmented plane wave (FP LAPW) code 
WIEN2k and the Vienna ab initio simulation package 
VASP 17 , using the projector augmented wave scheme^. 
In the FP LAPW calculations the Fe/Ir(001) systems 
were modeled by the eleven-layer repeated slab with ss 
10 A vacuum in between. All slabs were symmetric with 
respect to the middle layer. We allowed the relaxations 
for top three layers, the remaining interlayer distances 
were fixed to the bulk values of 1.92 A. For VASP re- 
peated asymmetric slabs with seven layers Ir and a single 
Fe mono-layer on one side and also symmetric slabs with 
eleven substrate layers and Fe mono-layers on both sides 
were used, which were separated by at least 19 A vacuum. 
All layer distances have been relaxed, and turned out to 
be essentially the same for both setups. Two DFT poten- 
tial approximations have been employed: the local den- 
sity approximation^ 5 - (LDA) and the generalized gradient 
approximation^ (GGA) for WIEN2k; the GGA accord- 
ing to Perdew and Wang (PW91)2 7 - as well as the LDA as 
given by Perdew-Zunger-22. (Ceperly-Alder) 2 ^ for VASP. 
We have tested NM-, FM-, and c(2x2)-AFM magnetic 
arrangements, all performed in the c(2x2)-structure for 
a reliable comparison of total energies. Technically, we 
have used a Brillouin zone sampling with 21-36 special k- 
points in the irreducible two-dimensional wedge. The dif- 
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ference between input and output charge density in the 
final iteration was better than 0.1 me a.u.~ 3 . The to- 
tal force on single atoms was in every case smaller than 
lmRy/bohr. 

In all TB-LMTO-SGF calculations the LDA approxi- 
mation and experimental layer relaxations^ were used. 
The vacuum above the overlayer was simulated as usual 
by empty spheres (ES). Electronic relaxations were al- 
lowed in four empty spheres adjoining the overlayer, the 
overlayer itself, and in five adjoining Ir substrate layers. 
This finite system was sandwiched selfconsistently be- 
tween an frozen semiinfinite fcc-Ir(OOl) bulk and the ES 
vacuum-space including the dipole surface barrier. Ad- 
ditionally to NM-, FM- and c(2x2)-AFM configurations, 
we have studied DLM-arrangements. 

In the framework of the TB-LMTO-SGF method the 
exchange integrals jf°' F ° between sites i,j in the mag- 
netic overlayer may be expressed as follow a 15 ' 24 
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Here, the trace extends over s— ,p—, d—, /-basis set, the 
quantities Af c are proportional to the calculated ex- 
change splittings, and the Green function • describes 
the propagation of electrons of a given spin (a =1,1) 
between sites i,j. It should be noted that both the di- 
rect propagation of electrons in the magnetic overlayer 
and the indirect one via the Ir-substrate are included in 
Eq. (JTJ on an equal footing. Finally, the energy inte- 
gration extends over all occupied valence states up the 
Fermi energy Ep which is technically performed by inte- 
grating over the contour C in the complex energy plane. 
For more details see Ref. [l5|. Once the exchange inter- 
actions were known, we constructed a two-dimensional 
(2D) classical Heisenberg Hamiltonian to describe the 
magnetic behavior of the Fe-overlayer on a non-magnetic 
fcc-Ir(OOl) substrate 



Fc,Fc 
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In Eq. ([5]) , denotes the orientation of the Fe-magnetic 
moment at the site i. By construction, the value of the 
corresponding magnetic moment is included in the def- 
inition of jfj e ' F °, and positive (negative) values denote 
FM (AFM) couplings. An early study of exchange inter- 
actions in fcc-Fe,Co/Cu(001) systems is found in Ref. H, 
whereas some recent estimates of exchange integrals for 
a bcc-Fe/W(001) overlayer were obtained either by a su- 
percell approac h 10 ' 30 or by an approach closely related to 
the present on o 31 ' 32 . 

A Green-function approach like the present one to cal- 
culate exchange interactions has a particular advantage 
over a supercell approac h 10 ' 30 : Exchange interactions can 
be evaluated easily and reliably even for disordered over- 
layers and partial coverages. 



III. RESULTS AND DISCUSSION 

A. Structure and magnetism 

The results of the total energy calculations are sum- 
marized in Tables I- IV. 

TABLE I: Calculated ( LDA/GGA, WIEN and VASP, respec- 
tively) and experimental (LEED) interlayer distances dij 
between top three sample layers (1-Fe overlayer, 2-top Ir layer, 
3-second Ir layer) for fcc-Fe/Ir(001) in the nonmagnetic (NM), 
ferromagnetic (FM), and c(2x2)-antiferromagnetic (AFM) 
states. Ht,-AFM displays the influence of 0.5 ML hydrogen 
adsorbed on favourable bridge positions. The interlayer (001)- 
distance in the bulk iridium is 1.92 A. 





di2 [A] 


d 2S [A] 


d si [A] 


LDA 


GGA 


LDA 


GGA 


LDA 


GGA 


NM 


WIEN 
VASP 


1.52 
1.51 


1.61 
1.58 


1.95 
1.99 


2.00 
2.05 


1.86 
1.88 


1.93 
1.94 


FM 


WIEN 
VASP 


1.64 
1.60 


1.78 
1.76 


1.91 
1.94 


1.95 
1.98 


1.88 
1.89 


1.94 
1.97 


AFM 


WIEN 
VASP 


1.59 
1.55 


1.69 
1.66 


1.93 
1.97 


1.98 
2.02 


1.88 
1.88 


1.93 
1.94 


H 6 -AFM 


VASP 


1.58 


1.67 


1.97 


2.01 


1.88 


1.94 


LEED 




1.69 


1.96 


1.91 



In Table I we present results of the structural mini- 
mization and compare our results with experiment. We 
found that the ground state is antiferromagnetic and that 
the interlayer distances obtained for this order in GGA 
approximation agree very well with the results of LEED 
structure analysis^. For all calculations we used the ex- 
perimental lattice constant a=3.84 A in layers, which lies 
between the calculated LDA bulk value (3.81 A) below 
and the GGA (3.87 A) above. For this reason, the calcu- 
lated substrate interlayer distances are slightly underes- 
timated in LDA (ss 0.04 A) and overestimated in GGA 
(~ 0.03 A), as the system tries to keep its respective 
equilibrium volume. The magnetic state influences the 
equilibrium surface geometry considerably at the Fe/Ir 
interface, while the changes in the substrate are less pro- 
nounced. A possible contamination with hydrogen, due 
to the preparation process should be hard to detect, as 
the hydrogen induced changes in the spacings are rather 
small and around the experimental error limit. The dif- 
ferences for first interlayer distance between the FM and 
c(2x2)-AFM configurations amount to ~ 0.10 A, while 
the NM Fe/Ir spacing is even smaller by 0.08 A. The cal- 
culated top layer relaxations (about 12%) is smaller than 
the one obtained for the similar bcc-Fe/W(001) system 
(about 14-19 

The first Ir-Ir distance is slightly expanded with re- 
spect to its bulk value, where one has to keep in mind 
that the bulk spacing is enhanced for GGA and decreased 
for LDA. The next spacing is reduced leading to a os- 
cillatory pattern of interlayer distances found in many 
metallic systems. 

The calculated work functions are presented in Table 
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TABLE II: Calculated work functions 5> in eV for Ir(OOl) (ne- 
glecting a possible lateral reconstruction) and of Fe/Ir(001) 
in various magnetic states. Symbols Ir, NM, FM, AFM, DLM 
Hi,- AFM, and Hj,-FM denote respectively, the Ir(OOl) surface, 
and nonmagnetic, ferromagnetic, c(2x2)-antiferromagnetic, 
and disordered local moment states of the Fe/Ir(001) over- 
layer, as well as the FM and AFM states with 0.5 ML hydro- 
gen adsorbed on bridge positions. The values correspond to 
the respective calculated relaxed geometries (see Table |TJ) for 
WIEN/VASP, and to the experimental one±i for LMTO. The 
experimental value for the Ir(001) surface is 5.67 eV— . 



<E>[eV] 


Ir 


NM 


FM 


AFM 


DLM 


H(,-FM 


H 


b-AFM 


WIEN-GGA 


5.65 


4.86 


4.38 


4.45 










VASP-GGA 


5.62 


4.82 


4.29 


4.37 




4.71 




4.65 


WIEN-LDA 


5.92 


5.14 


4.58 


4.66 










VASP-LDA 


5.89 


5.13 


4.59 


4.66 




5.01 




5.04 


LMTO-LDA 


6.22 


5.05 


4.79 


4.76 


4.67 









TABLE III: Calculated stabilities (in mRy/Fe atom) of the 
various magnetic phases of Fe/Ir(001) as obtained by the 
WIEN, VASP, and LMTO codes for the unrelaxed geome- 
try (dp _i r =1.92 A). The ferromagnetic ground state has the 
lowest energy and serves as the point of reference. 

A[mRy] NMlda AFM lda DLM LDA NM G ga AFMgga 



WIEN 


49.1 


2.2 




62.7 


2.2 


VASP 


44.5 


3.8 




60.4 


3.9 


LMTO 


42.5 


5.1 


2.1 







II. The experimental value of 5.7 eV for the fcc-Ir(OOl) 
is reasonably reproduced by the present calculations (no 
surface reconstruction). Our calculations show that the 
iron overlayer reduces the sample work function by ~ 
1 eV (no experimental data have been found). Hydro- 
gen adsorbed on the overlayer should increase the work 
function by « 0.4 eV. The calculated values for the over- 
layers are close to the work function of bcc-Fe (about 
4.5 eV) obtained for a polycrystalline sample 33 . While 
the full potential VASP and WIEN codes agree very well 
with each other, the TB-LMTO-SGF approach slightly 
overestimates the values. 

The results of magnetic stability calculations are pre- 
sented in Tables III - V. Different theoretical approaches 
were used and compared in Tables for the unrelaxed ge- 
ometry (Table III) as well as for the realistic, relaxed 



TABLE IV: Calculated stabilities (in mRy/Fe atom) of the 
various magnetic phases of Fe/Ir(001) as obtained by the 
WIEN, VASP, and LMTO codes for relaxed geometries. The 
stabilities correspond to the respective calculated relaxed ge- 
ometries (see Table HJl for WIEN/VASP, and to the experi- 
mental oneM- for LMTO. The antiferromagnetic ground state 
has the lowest energy and serves as the point of reference. 

A[mRy] NMlda FMlda DLM lda NMqga FMqga 



WIEN 


25.8 


7.8 




38.4 


5.3 


VASP 


22.1 


8.6 




38.7 


5.5 


LMTO 


27.7 


5.0 


0.8 







TABLE V: Influence of 0.5 ML adsorbed hydrogen on the 
calculated stabilities (in mRy) of the various magnetic phases 
of Fe/Ir(001) as obtained by VASP for relaxed geometries (see 
Table [TJ . The antiferromagnetic ground state has the lowest 
energy and serves as the point of reference. 

~A[nrRy] NMlda FM lda NMgga FMgga" 

clean 22.1 8.6 38.7 5.5 
bridge-H 14.4 3.7 29.5 2.4 



case (Table IV) including the possibility of residual ad- 
sorbed hydrogen (Table V). It should be noted that one 
may only compare different LDA or GGA energies di- 
rectly to each to other. Our calculations clearly show, 
that the nonmagnetic case can be safely excluded. All 
models with local magnetic moments have a substan- 
tially lower total energy. The most striking result, ob- 
tained by all methods, is the fact that while the FM is 
the ground state for an unrelaxed geometry, the layer 
relaxations stabilize the c(2x2)-AFM phase. This is in 
a striking contrast to the closely related bcc-Fe/W(001) 
casein where the antiferromagnetism of bcc-Fe/W(001) 
is robust with respect to the structural relaxations. In- 
terestingly adsorbed hydrogen also does not change the 
picture as evident from Table V. Differences only get 
smaller, but the general trend is preserved. However, 
strictly speaking, the c(2x2)-AFM may not be the true 
ground state of fcc-Fe/Ir(001). Similar to bcc-Fe/W(001) 
one should consider other possibilities, e.g., the p(2xl)- 
AFM or even some non-collinear configurations, such as 
spin-spirals or a chiral stated. To shed some light on 
this issue we included in Tables III and IV the results 
of DLM calculations as performed in the framework of 
the TB-LMTO-SGF approach. The sufficient reliability 
of the TB-LMTO-SGF approach for the present purpose 
is confirmed by a comparison with accurate full potential 
calculations for NM-, FM-, and c(2x2)-AFM configura- 
tions in both the ideal and relaxed geometries. This also 
justifies the use of the TB-LMTO-SGF approach to ob- 
tain exchange interactions in the Fe-overlayer in the fol- 
lowing. It is quite obvious that, if the DLM state is the 
ground state compared to the NM- and FM-states, then 
a more complex, AFM-like state can exist and this con- 
clusion can be reached without performing calculations 
for many possible candidates. But clearly, this fact does 
not render the necessity of searching for the true ground 
state of the system obsolete, but rather represents a re- 
liable qualitative indicator for a more complex magnetic 
state of the system. 

One can speculate about the structural origin of such 
an AFM order. As already mentioned, both bcc(001) 
and fcc(001) have a common square-lattice structure of 
magnetic atoms which differ by the ratio d : a^, where 
<iL is the lateral lattice constant and d is the layer spac- 
ing. For bcc(001) we have d = Ol/2 = a/2 (where a is 
the bulk lattice constant) while for fcc(001) d = a/2 and 
cil = a /V2 leading to a larger d : a^. The reduction 
of the d : a,L ratio stabilizes the AFM-/DLM-state for 
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this fcc(OOl) surface. This ratio is sufficiently small for 
comparable bcc(00f ) surfaces even in the unrelaxed ge- 
ometry (see e.g. Fe/W(001) 8 ). ft is well-known that the 
exchange interaction between Mn spins becomes antifer- 
romagnetic for smaller distances, a trend we observe here 
for exchange interactions in the Fe-overlayer as detailed 
below. This is a first strong indication that indirect inter- 
actions of Fe-spins via the Ir-substrate play an essential 
role for the fcc-Fe/Ir(00f ) magnetism. A dominant char- 
acter of indirect interactions as compared to direct ones 
between Fe-spins in the overlayer will result in a strong 
dependence of exchange interactions on the Fe-Ir inter- 
layer distance as we shall see below. 

B. Densities of states and exchange interactions 
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FIG. 1: Total layer-resolved densities of states (DOS) for the 
fcc-Fe/Ir(001) overlayer and experimental layer relaxations 11 
based on the LMTO approach. In the case of the Fe-overlayer 
the total DOSs are additionally split into majority (dashed 
lines) and minority (dotted lines) contributions: (a) the FM 
case, (b) the DLM case, and (c) the c(2x2)-AFM case. Only 
one spin-orientation (denoted as Fe + ) is plotted for the DLM 
and AFM cases. Symbols Vac, Fe, ST, S-2, and bulk denote 
the first vacuum layer, Fe-overlayer, first and second substrate 
Ir-layers, and the fee bulk Ir host DOSs, respectively. The 
Fermi energy is shifted to the energy zero. 

In Fig[T]we present the layer-resolved densities of states 
(DOS) for the various magnetic states. Compared to the 
Ir bulk, the most important feature observed for all cases 
is an extra contribution to the overlayer DOS around the 
Fermi energy due to the minority Fe-states. The large 
exchange splitting of majority and minority Fe-states is 
due to the enhanced overlayer magnetic moment (about 
2.65 [is in ah cases) which illustrates the rigidity of the 



Fe-moment with respect to changing spin-orientations. 
The large Fe moment is due to the large lateral overlayer 
lattice constant as given by the Ir-substrate as well as 
to the reduction of coordination number at the surface 
typical for overlayer systems. There is also a relevant 
extra peak in the DOS in the vacuum close to the sam- 
ple surface which gives a possibility to detect it in STM 
measurements. We also observe a strong reduction of the 
Fe-overlayer imprint on the deeper Ir-substrate layers: 
already the second Ir-substrate layer is almost bulk like. 
There is a small induced moment on the first-substrate 
Ir-layer of the order of 0.1 [Lb while other induced mo- 
ments are strongly damped in an oscillatory manner (the 
Friedel-like oscillations) into the Ir-substrates and their 
values are of the order of 0.01 \ib and smaller. It should 
be noted that induced moments in the Ir-substrate and in 
the vacuum are much smaller and more strongly damped 
in both substrate and vacuum for the c(2x2)-AFM case 
as compared to the FM case, while they even collapse 
to zero for the DLM state. The induced moment in the 
substrate/ vacuum are more than two-times smaller as 
those in the bcc-Fe/W(001)^. The above results are 
confirmed by full-potential slab-model calculations using 
both WIEN and VASP codes. 

Exchange interactions J Fe,Fe (d) have been determined 
by the TB-LMTO-SGF method for Fe/Ir(001) as a func- 
tion of the interatomic Fe-Fe distance d for both un- 
relaxed and relaxed geometries, and FM-, DLM-, and 
AFM-reference states. Additionally we present results in 
the DLM state for a simple layer-relaxation model where 
only the Fe-Ir interlayer distance is reduced from the un- 
relaxed value (1.92 A) to values of 1.82 A, 1.72 A, and 
1.62 A corresponding to a reduction of 5 %, 10.5 %, and 
16 %. These results are shown in Fig. O 

Calculated exchange integrals are rather similar, e.g., 
all leading to the AFM interactions for the relaxed case 
irrespective of the reference state. The problem of the 
choice of reference state for estimate of exchange inte- 
grals in the present context was also addressed in Ref . [34]. 
We have chosen the DLM reference state (with self- 
consistently calculated spin-polarized potentials and with 
corresponding modification of Eq. |[T])) for the magnetic 
stability study below as it assumes no magnetic ordering 
and for our purposes is thus most suitable. We admit that 
to estimate e.g. the critical temperatures other choices 
may be more suitable. 

With increasing layer relaxations we observe a clear 
tendency towards dominating AFM interactions which 
stabilize the AFM-like state in the overlayer. The domi- 
nating role of indirect interactions between Fe-atoms via 
the Ir-substrate is obvious: the only varying quantity is 
the Fe-Ir distance and thus the Fe-Ir hybridization. Re- 
sults for experimental layer relaxations (inset in Fig. [5]) 
are between model cases of 1.72 A, and 1.62 A. The fact 
that strong AFM coupling in the layer-relaxed case were 
obtained from a reference FM state indicates the robust- 
ness of the AFM order (more precisely, of a more complex 
magnetic state) for the relaxed Fe/Ir(001) overlayer. 
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FIG. 2: Exchange interactions among Fe-atoms in the 
Fe/Ir(001) overlayer for the unrelaxed case and model layer 
relaxations evaluated as a function of the reduced interatomic 
distance (d/a), where a denotes the lattice constant. Num- 
bers attached to symbols indicate the reduction of the Fe-Ir 
interlayer distances in % as compared to the bulk value of 
1.92 A. The inset shows exchange interactions for AFM and 
DLM state for the experimental layer relaxations^. All re- 
sults were obtained assuming the DLM-reference state. 
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FIG. 3: Lattice Fourier transformation of the real-space ex- 
change interactions J Fc,Fc , J(q||), for the ideal unrelaxed ge- 
ometry, the experimental layer relaxations" as well as for 
three model Fe-Ir layer relaxations (same as in Fig. [2]) which 
were obtained for the reference DLM state of the Fe/Ir(001) 
overlayer .Here, q||=X=27r/a £ (l/2, 1/2), q||=f=27r/ai,(0, 0), 
and q||=M=27r/o£(l,0). 



To further investigate this point we present in Fig. [3] 
the result of the lattice Fourier transformation of ex- 
change integrals {1} as obtained for the DLM reference 
state. It should be noted that the DLM state is neither 
the ground state for the ideal geometry nor for the exper- 
imental layer-relaxed model. On the other hand, possible 
magnetic phases of the 2D-Heisenberg model ([2]) can be 
obtained by studying its stability with respect to the pe- 
riodic excitations. A similar approach was successfully 
used in the study of the complex magnetic stability of 
bcc-Eu: starting from the FM reference state, a proper 
spin-spiral ground state was obtained in good agreement 
with the experiment^. Due to the sign convention in 
eqaution |2]) the maximum of J(q||) corresponds to the 
ground state (the energy minimum). It is obvious that 
the ground state for the unrelaxed model is the ferromag- 
netic state (the maximum of J(q||) is obtained for q|| = 
2ir/ciL(0, 0). With reduced Fe-Ir interlayer distance we 
observe a quick decrease of the stability of the FM state 
and the ground state is found for the ordering vector 
q||=27r/a L (l/2, 1/2) which corresponds to the c(2x2)- 
AFM state (for the Fe-Ir interlayer distance d=1.82 A or 
reduced by 5%). If the Fe-Ir interlayer distance further 
decreases (by 10.5% and 16%) the stability of the FM 
state further decreases and a new, complex spin-spiral 
like magnetic state becomes more stable as compared 
to the c(2x2)-AFM state (ordering vector on the line 
X-M in the irreducible surface Brillouin zone). On the 
other hand, the c(2x2)-AFM state becomes more and 
more stable as compared to the FM state in accordance 
with the total energy calculations. Also, both for the 
experimental layer-relaxation modelii and for the model 
with the largest reduction of the Fe-Ir interlayer distance 
(by 16%) the p(2xl)-AFM state (the ordering vector 
q||=27r/a_L(l, 0)) has lower energy as compared to the 
FM- and c(2x2)-AFM states. 

In order to advance the understanding of the present 
substrate induced coupling a comparison with an other- 
wise similar 4d substrate would be beneficial. Rhodium 
crystallizes also in the fee structure, has the same number 
of valence electrons (9) as iridium and a similar lattice 
constant (3.80 A), while the value of the work function, 
5.11 eV— , is smaller as well as the spatial extent of the 
4d-wave functions of Rh as compared to 5d-wave func- 
tions of Ir. The results of a similar study of the mag- 
netic stability for fcc-Fe/Rh(001) overlayer as a function 
of the Fe-Rh interlayer distance will be presented in the 
following. The same relative interlayer reductions as for 
fcc-Fe/Ir(001), namely by 5 %, 10.5 %, and 16 % will be 
used. 

The results are summarized in Fig. [4] The FM state 
is again the ground state for the unrelaxed case, and we 
observe a similar Fe-Rh interlayer distance reduction ef- 
fect on the magnetic stability as for fcc-Fe/Ir(001). The 
reduction of the Fe-Rh interlayer distance by about 5 % 
seems to be the point where the FM state is no longer the 
ground state and the c(2x2)-AFM state is stabilized (the 
ordering vector =27r/aL(l/2, 1/2)). However, with in- 
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FIG. 4: Lattice Fourier transformation of real-space ex- 
change interactions Jff' Fc , J(q||), for the unrelaxed geome- 
try, and three model Fe-Rh interlayer relaxations for the fcc- 
Fe/Rh(001) system. The identical relative interlayer distances 
were used for both the fcc-Fe/Ir(001) system and for the DLM 
reference state. 



creasing reduction of the Fe-Rh distance the more com- 
plex AFM state is stabilized until for an interlayer re- 
duction of 16 % this state become the ground state (see 
an indication of this case in the Fig. 4: maximum in the 
neighborhood of the q||=27r/ai(l, 0) ordering vector rep- 
resenting the p(2x 1)-AFM state). Again the AFM state 
is stabilized but for a relatively larger reduction of the 
Fe-Rh interlayer distance as compared to fcc-Fe/Ir(001), 
but otherwise there is a large similarity between the two 
over layer systems. The larger reduction is in agreement 
with previous total energy study-i^. 



IV. CONCLUSIONS 

We have investigated the experimentally prepared 
Fe/Ir(001) system by a combination of different first- 
principles methods for both supercell slab geometries 
(WIEN and VASP codes) and semi-infinite boundary 
conditions (TB-LMTO-SGF codes). Using the latter ap- 
proach we investigated the magnetic phase stability of 
the system by calculating the exchange interactions be- 
tween Fe-atoms in the overlayer. The following conclu- 
sions can be drawn: (i) Calculated relaxed geometries 
agree well with available experimental data. A better 
agreement is obtained for the Fe-Ir distance using GGA 
rather than LDA potentials. The calculated geometries 
depend on the magnetic state (non-magnetic, FM, and 
AFM states) but the most important result is the re- 
duction of the Fe-Ir interlayer distance as compared to 
the bulk value by about 12 % for the AFM order. A 



possible residual H contamination on the overlayer has 
practically no effect on the distances; (ii) The local Fe- 
moment is enhanced to about 2.65 fj,B as compared to 
its canonical value of 2.15 [Ib in the bcc Fe-metal. The 
enhancement is due to both the enlarged lateral lattice 
constant of the Fe-overlayer on fcc-Ir(OOl) and the re- 
duction of nearest-neighbors there. The work function of 
the system with Fe-overlayer is reduced more than 1 eV 
as compared to the value found for pure Ir-surface. Hy- 
drogen on the overlayer increases the work function by 
« 0.4 eV; (iii) Both full potential supercell methods found 
the c(2x2)-AFM state to be stable as compared to the 
non-magnetic and ferromagnetic states for the correctly 
relaxed structure. On the other hand, they find an FM 
ground state for a unrelaxed geometry. Hence, the sta- 
bility of the AFM state is induced by the substrate via 
layer-relaxations. The above results were confirmed both 
qualitatively and quantitatively by the TB-LMTO-SGF 
approach assuming a scmiinfinitc sample geometry and 
experimental layer relaxations as obtained in LEED; (iv) 
The related fcc-Fe/Rh(001) system behaves similarly to 
fcc-Fe/Ir(001), although the tendency towards the forma- 
tion of the AFM state upon reduction of the Fe-Rh in- 
terlayer distance is weaker. We found an indication that 
a p(2x 1)-AFM state is stabilized for larger reductions of 
Fe-Rh interlayer distances; (v) A detailed study of the 
magnetic stability based on the 2D-Heisenberg Hamilto- 
nian derived from the first-principles total energies con- 
firms the stability of the c(2x2)-AFM state as compared 
to the FM for relaxed model but indicates that more com- 
plex, spin-spiral like state can be stabilized by reducing 
the Fe-Ir interlayer distance. Since Ir is a heavy element, 
the spin-orbit interaction can have non-negligible effect 
on the exchange interactions (in analogy to the Mn/W 
casei), and can lead to chiral magnetic order induced 
by Dzyaloshinskii-Moriya interaction; and (vi) We have 
shown that increasing Fe-Ir hybridization stabilizes the 
AFM state. This seems to be a robust effect to war- 
rant similar stabilization on a much more complex re- 
constructed Ir(001) surface. However, such a study goes 
beyond the subject of the present paper. 

Experimentally no-magnetization was found for iron 
overlayers thinner than ss 4 monolayers on fcc-Ir(OOl) 11 
and s» 6 monolayers on fcc-Rh(OOl) — . Our study indi- 
cates that this is not the result due to disappearing of the 
local iron magnetic moments in the top surface layers but 
rather due to a more complex " antiferromagnetic-like" 
order in the fee- iron overlayer. 
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